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In a recent work, factorization breaking interactions and T-odd distributions have been analyzed 
from the point of view of statistical mechanics. One of the central points was the process leading 
from the coherent initial hadron state to an incoherent set of short-range parton states, in presence 
of factorization breaking interactions. This time-evolution was named "partonization" , or also 
"decorrelation" , to distinguish it from an ordinary hadron-quark-spectator splitting vertex. Here an 
example is shown for a correlation amplitude describing an extreme case for such a process, where 
the initial state phases are lost, but the probability distributions are perfectly conserved. This 
requires introducing degrees of freedom that are not normally considered in a factorized scheme, 
Oh' and that must disappear at some stage of the calculations for the factorization scheme to survive. 

CD ' The explicit time-evolution character of the decorrelation process disappears when these degrees of 

m : 

freedom are integrated, but the integrated correlator presents modified analytical properties as a 
consequence. 
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I. INTRODUCTION 
A. General 



The study of the so-called T-odd distribution functions has made it necessary to reconsider some of the basic 
properties of the parton model[lJ and of the related factorization schemes (for the points concerning us here, see 
refs.@, Q). We may roughly identify two main research lines on T-odd distributions. 

One starts with the first T-odd "official" distribution, the Sivers function^] in 1990, followed by the Boer-Mulders- 
Tangerman function 0, 0, 0]. The former was used to explain single spin asymmetries [8], the latter to explain unpo- 
larized Drell-Yan azimuthal asymmetries H|. 

The other research line is olde r 13, [lOL and is centered on the study of twist-3 (with respect to qr) effects. It has 



been very recently demonstrated [12], |13|], that the two things (qr— twist-3 contributions and Sivers function) describe 
the same phenomenon in different gy— regions, and produce effective distributions with the same Q m — dependence in 
Drell-Yan (the latter describes objects that are twist-3 with respect to qx, but in Drell-Yan qx and Q are independent). 

For completeness, we should note that singl e sp in asymmetries and unpolarized azimuthal Drell-Yan asymmetry 
have also been interpreted via soft mechanisms [lj, although still according with schemes that may be (qualitatively) 
C3 ' considered "T-odd" for the role played by initial state interactions. 

In nuclear physics a time reversal odd structure function, the so-called "fifth structure function" , was introduced 
m 1983 by Donnelly[Il and much later modeled [l6l [ItJ to describe normal asymmetries in A(e, e'p) quasi-elastic 
scattering (for reviews see [18jV Although the formalism presents differences, several connections with the partonic 
T-odd functions are present. A point that is exploited later is the observation [l7j that a key passage for producing a 
nonzero fifth structure function is giving absorption features to the rescattering terms. 

In hadronic physics, a key role is played by the gauge restoring field operator entering the definition of parton 
distribution [2j ■ In an approach where factorization is assumed from the very beginning, this field is neglected, but 
in such case leading twist T-odd distribution functions are forbidden (l9| by general invariance principles. After a 
subtle analysis [13] had shown that proper taking into account the gauge factor implies leading twist corrections in 
small- a; deep inelastic scattering, an explicit mechanism was shown to produce a nonzero T-odd distribution in a QCD 
framework [2 ll l22j]. In [23j the existence problem was systematically fixed by analyzing the time-reversal properties 
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of the gauge factor, and in [2J, |25[ the pitfalls associated with the limiting behavior of the gauge fields have been 
examined in depth. Serious efforts has been undertaken [l^, [H, HH to rewrite factorization rules for processes where 
partonic distributions depend on (soft values of) transverse momentum and the gauge factor is assumed to play a 
role. In general however, the impossibility to neglect the gauge fi eld obliges one to reconsider the correctness, or at 
least the effectiveness, of the assumed factorization [27l [28L T29LT30I] . 

Several authors have given contributions to the huge development of this field in the last few years. Models (3ll. 
HH HH H3 and studies (35l. l36l. l37l. l38l l39l. l4fl l4lj around T-odd distribution functions have been produced. Recently 
phenomenological parameterizations of the Sivers function for qua rks [42 , |43L |44|. |45|| and gluonsTfj have been deduced 
from available data[U IH, HjJ and also parameterizations 0, l5l| of the Boer-Mulders-Tangerman function have 
been produced. Quite a few experiments aimed at the measurement of T-odd functions are planned for the next ten 
years ll^Ill [Hi, i|. 



B. Starting hypothesis and aims of this work 

In a previous work a picture was presented of T-odd functions as a result of a statistical process associated with 
initial state interactions between two hadrons. This transformed the hadron structure, that initially was characterized 
by long-range correlations ("hadronic phase") into an incoherent collection of short-range objects describing a quark 
in the quark-parton model ("partonic phase"). So, the presence of a nonzero T-odd function was related to the decay 
properties of the initial state wavcfunction, and to a singular rate of increase of the entropy of the system. 

A point that was not discussed in that work was the way this decorrelation process could show information about 
itself in the two-point amplitude describing a semi-inclusive measurement, in particular a Drell-Yan measurement (the 
so-called correlator of the process, see next section). The presence of such information in the correlator is essential to 
have a connection between the physical phenomenon (the decorrelation process) and the possibility to observe it in 
the distribution functions. 

The main section of this work is section V, where we show an example of a correlation amplitude describing a 
partonization/decorrelation process. This two-point amplitude contains in itself an initial pure state, a final incoherent 
probability distribution, and the transition between the two regimes. 

A peculiarity of the decorrelation process is brought to its extreme by the chosen correlator: the hadron-to-parton 
transition is "ideal" , i.e. although phases are destroyed in the transition, probabilities are perfectly conserved. We 
will speak of "ideal decorrelation/partonization process". 

The description of a process involving factorization-breaking interactions (causing the decorrelation) requires in- 
troducing extra degrees of freedom in addition to those associated with x and fey, the quark longitudinal fraction 
and transverse momentum. There is nothing exotic in this. Any model considering rescattering introduces such extra 
coordinates, and several of the above quoted references do this. In the following we will name "extra" degrees of 
freedom those that are not present in a standard factorization regime. 

Let £ = P+Z- be the spacetime coordinate conjugated with the quark longitudinal fraction x (see section II for all 
the definitions), and let G(£) be the correlation amplitude in a factorization-respecting approach. Let r\ be an extra 
coordinate introduced by rescattering. A G(£) correlator that at least formally respects factorization must originate 
in an integration process of the form G(£) = J dr)G'(^,Tf). 

We will consider a decorrelation process that is evident in the pre-integration amplitude G'(£, 77), but not in G(£). 
We will show that G(£) still keeps in itself a track of the decorrelation process, that a hard probe is able to detect 
in the form of a T-odd additional distribution function. This track has to do with the analytical properties of the 
correlator. 

Distribution functions are calculated from the imaginary part of the correlator. This imaginary part derives from a 
real particle production cut, and the cut position has an infinitesimal imaginary shift iO from the real axis associated 
with causality. We will show that the imaginary shift becomes finite: iO — ► in, n > 0, when the proposed correlator 
is averaged over extra degrees of freedom. 



In another work 57j it is shown that such a finite shift may imply T-odd observable corrections to the distribution 
functions. Here we only discuss the formation of the finite imaginary shift. 

We stress that in the following nothing is demonstrated about factorization. But we require the assumption 
that factorization breaking is a small effect, because this allows for simplifying approximations. In concrete terms, 
"small effect" means that the parameter n (appearing in the following associated with rescattering corrections to the 
correlator) is small. 
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C. Additional degrees of freedom 

The relevant distribution functions are extracted from a two-point correlation function G(l,2). This describes the 
propagation of a quark hole created in a hadron in the spacetime position "1" to the spacetime position "2". 

In presence of factorization, G depends on "1" and "2" via the difference "1-2" only (the absence of rescattering 
implies translation invariance). In addition, one of the two light-cone projections of "1-2" must play no role. 

There are several degrees of freedom that one may imagine to acquire a role when factorization breaking interactions 
are included. First, the ones related with "1+2", or with the neglected light-cone component, or with the position 
and internal structure of the target hadron. 

To simplify things, we limit here to one extra degree of freedom. This is sufficient to understand how decorrelation 
develops and how it leads to dissipative effects when averaged. We neglect z+ and the other possible variables, but 
include "1+2" (specifically: z-(l) + z_(2)) in our calculations. 

II. GENERAL DEFINITIONS AND NOTATIONS 

We consider a Drell-Yan process taking place between two colliding hadrons with momenta P M and P^, in their 
center of mass frame. The parton we will examine in the following is a quark belonging to the former hadron. 

To select leading twist terms, we will use the infinite momentum limit P+ — ► oo. In this limit (P + P') 2 ~ P+P'_ 
~ P + 2 is infinite, and for finite x and x also Q 2 = xx(P + P') 2 is infinite. 

The relevant amplitude considered in this work is the traditional one where a hole is created in a hadron state by 
extracting a quark/ antiquark in a point z'^, the "hadron+hole" set propagates across the cut through real states only, 
and in z"^ the hole is filled and the original hadron state is restored by reinserting the quark in its place. All the 
observable quantities interesting us are extracted from the imaginary part of this amplitude, that in this paper we 
simple name (7(1, 2). Where necessary, the spacetime coordinates corresponding to "1" and "2" are explicitly written. 
Variables corresponding to "1" are indicated as z', k' etc, variables corresponding to "2" with z", k" etc. 

The extracted quark has momentum k^ with lightcone fc+ = xP+ and negligible fc_. The impact parameter b is 
the space vector conjugate to kx- The spacetime displacement of the hole from z'^ to z 1 ' is z^ = z'^ — z 1 ' = (z_,z+,b) 
w (z-,0,b). The fourth coordinate z + plays no role and is not explicitly reported in the following. 

Any relevant distribution function q(x) is the Fourier transform with respect to z_ and & of a correlation function 
g{z^,b). Leading twist effects are naturally selected by writing the Fourier transform with respect to the scaled 
variables x and £: 

i = p+z-, (i) 

q(x,k T ) = P+ J e- txP + z -e tkTb g(z-,b)dz-db = J e - lx ^e lkTb G(Cb)d(db. (2) 

G(£, b), £, x, etc will be named "scaled" quantities. The range of useful values of £ remains finite ~ 1/x when P+ 
— > oo. In this limit the function (7(£) does not need to be singular in the origin to produce a nonzero q(x). So, in the 
following the discussion is in terms of the scaled quantities. 

G is here used as a function and not a matrix, so it must be read as a projection of the spinor correlation matrix Gij 
over a suitable operator, e.g. 7 + . The 7 + projection selects simultaneously the main unpolarized quark distribution 
and the Sivers T-odd function. 

The two functions qrx, kx) and G(£, b) are related by double Fourier transformation with respect to x/£ and kx/b. 
In the following, we will mostly work on the mixed representation function /(x, b): 

f(x,b) ee / e-^G^b)^. (3) 



G and / depend on longitudinal and transverse variables. 

Where not necessary to resolve ambiguity, the transverse variables will not be explicitly written, since the amount 
of nonstandard variables and constants n, A, etc needed in the following is rather large. So, e.g. the above eq|3]may 
simply assume the form 
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f(x) = / e-^G(Od( (4) 



where the presence of b is implicitly assumed on both sides. This should cause no ambiguity, since in this paper 
kx— integrated quantities do not appear: all the considered functions depend on both longitudinal and transverse 
variables. 

In the following, we will have to extract G(£) = G(£' — £") from the integration of a two-point amplitude G(£', £") 
over £' + £". The latter can be written in terms of two independent variables x' and x" via a double Fourier transform 
exp[— i(x'd' — x" £_")]. To separate degrees of freedom, we use the definitions 

e = z'-e, v = ^4^- (5) 



X = A = (6) 

then we have 

exp[-i{x'? -x"£')] = ex P [-i{X£ + At])} (7) 

III. LIGHT-CONE TRANSLATION AVERAGING OPERATOR 

In presence of a violation of translation invariance, we could define f A (x) as an average over 77 = (£' + £")/2: 



(8) 



and a distribution /(a:) could be defined as f A (x), if A is large enough to wash away the effects of 77— irregularities. 
To simplify notations we define 1 



[I A G}(0 = G(?,adv (9) 

assuming that the dependence on rj is completely lost when the average covers a large enough range A. We may 
rewrite the above eqEl 

f A (x) = IaJ e-**tG(Z',Z")dZ = J dCe^[l A G}(0- (10) 

Since the ^' + £" and — £" integrations can always be exchanged as we have just done, a translation invariant object 
may always be defined, for a large enough A. 

The factor 1/A is not necessary but guarantees a correct normalization of final results, according with later equations 
in this work. One may equivalently define I A as a sum operator, rather than an average operator. However, either 
in I A , or in the definition of the wavefunctions entering the correlation function, a normalization must be present. 
Formally, I A may be applied, redundantly, also in absence of factorization breaking effects. Then, the normalization 
1/A leaves the results unchanged, without the need of introducing further normalization factors in the wavefunctions. 

Instead of I A , we could use another averaging operator: 



1 The integration range may be also taken as symmetric from —A to 
have £' > £" due to intermediate real state propagation from £" to 
contained in this work are even with respect to r], one may limit to 



+A (from —00 to +00 below in the definition of Km)- Indeed, we 
but this does not touch r\ = (£' +£")/2. Since all the relations 
positive r\ only. 
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M r+oc 

ds 

-M Jo 



[K M G}(0 = 1 / d, / '"'•(7(i'. i "i (ID 



where s is the longitudinal fraction conjugated with rj. This may look more appropriate than I a, since it introduces 
the double Fourier transform eqEJ and the sum is over a momentum-like loop variable. However, Km is equivalent to 
I A, for N — A and M — ir/A. Exchanging the i] and s integrations, and applying the operator to G(rj, £) = G(£,£"), 

r°° r M r°° sin(Mn) r /M 

/ drjG^O / dse-^ = 2 / (ftygfofl K ' « 2 / £)• (12) 

Jo J-Af Jo *7 Jo 

So from now on we will use directly I a ■ 

IV. ENLARGED DELTA-FUNCTIONS 

If the problem is rj— independent, a S(x' — x") function associated with the conservation of the light-cone momentum 
is present (according with eqs. [6] and [7J a;' — x" is the variable conjugated with rj). In presence of interactions 
that violate x— conservation with \x' — x"\ ~ n, we need substituting S(x' — x") with a function that takes this 
nonconscrvation into account. 



Exploiting a technique developed in 17[ and [581 ] to roughly but quickly approximate the resummed effect of 
rescattering, we define an "enlarged 5— function" S n (x' — x"): 

(i) x' w x" within gaussian uncertainty |n|, 

(ii) Jq dx"8 n {x' — x") 1 (with exact equality for an infinite integration range). 

Such functions are normally used to define the ordinary S(x' — x") function via their n — » limit. 
Several equivalent representations for 5 n (x' — x") are possible, and in the following we will use different ones 
interchangeably, accordingly with mathematical convenience. Also, factors like \/3 or tt will be normally neglected. 
A useful representation is 2 

S n (A) = 



tt n 2 + A 2 

sK , w + K _ Ws ± T L Tn -i- T L rn (i3) 

Causality implies the substitution 

S n (A) -^S n+ (A). (14) 

that corresponds, in its Fourier transform, to 

exp{-n\£\) -> 9(£)exp(-nt;) (15) 

Since £ > in the relations interesting us, in the following S n must be meant as <S n +, i-e. one pole only is considered 
in the complex plane integrations. At the condition that n is small enough to justify the idea that the (real-axis) 
integral is dominated by the peak at A = 0, we may approximate the integral with the contribution from the pole x' 
= x + in of the function i5 n + (see (57j for a more precise justification of this): 

1 df 
dx'f(x')5 n (x-x') sa f(x + in) sa f(x) +in+ + 0(n 2 ). (16) 

dx 



From now on we use A = x' — x" , X = (x' + x")/2, see eq|6] 
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In the following, wc need the relation 3 

I A [5 lh m « S 1/A (A) (17) 
To show its validity, we use the another representation for an enlarged S— function: 

i r+\v\h 

<5 7/ ,(A) = - / e^dy (18) 
2lT J-\v\h 

When I a is applied to it, we get 

27r&(A) = iA[5 l/v (A)} = 1 f A d V f +Mh e">' A d V ' (19) 

A Jo J-Wh 

where we have defined 

n = j/A. (20) 

That S n is really an enlarged— 5 function may be demonstrated by direct integration of the previous equation: 

1 sin 2 (A/n) 
2™ (A/n) 2 



The right hand side function has n— independent A— integral = 1 and limiting behavior (for small n) like a delta 
function (its value is oo for A = and for any other A). So it is a fair representation for an n-enlarged <5 with 
precise width y/3n (from the small- A power series for it). 

V. A DECORRELATING CONFIGURATION 

We assume that initial state interactions, and consequently the correlator, depend on both £ and rj. We introduce 
a parameter 7 so that the effect of rescattering in terms of x— nonlocality decreases at increasing values of 77, to loss 
any effectiveness for 77 >> 7. This makes sense, since a large r\ means large values of £" or both (their difference 
£ may remain small) meaning that the target hadron is far and interaction is suppressed. 

A function of x that depends parametrically on 77 can be used to model an rj— evolution towards decorrelation. 

Let us use the limiting and normalization properties of 6 1 / v (x' — x")\ 

5 7 / ri (x — x') w S(x — x'), for 77 >> 7; (22) 
5 7 / v (x — x') w const, for \x — x' | << 7/77, (23) 

J dx'5 lh {x- x') Rj 1. (24) 



Wc define 



j dxip(x) 



tf(f) = / dxiP{x)e ixii , (25) 



3 The use of the sign "as" instead of "=" is due to (i) a constant factor of magnitude unity, (ii) two different representations are used for 
the functions 5 n appearing in the left and right hand sides of the equation, and the corresponding functions arc equal (in distribution 
sense) only in the n — > limit. 
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G(£',£") ee G(r),0 ee J dx'dx"^{x')^*{x")e^ x ' +x " ) ^'^" )/2 5 lh (x-x'). (26) 

dXdA + A/2)tp*(X - A/2)e iXS <5 7/j7 (A), (27) 
G has the following limiting properties: 

G(£',£") « G (0, /ortj « 7 ; (28) 

« Goo(0, for V » 7. (29) 

where 

G (0 ee *(e/2)**(-$/2), (30) 

Goo(f) s | dX|^(X)|V*«. (31) 

The former object is the correlator of a pure state, while the latter is the correlator of an incoherent overlap of 
plane waves exp(— ix£). Remarkably, in both regimes the probability for extracting a given plane wave exp(ix^) from 
the system is the same: |^(x)| 2 . 

So we get an "ideal partonization process" with respect to x. The shrinking width-parameter of the enlarged delta 
handles the transformation from the density matrix of a pure state to the one of a probabilistic distribution. 

The separation parameter of the two regimes is 7. For rj much smaller or larger than 7 we are in one of the two 
regimes. As far as G is not averaged with respect to 77, it contains explicitly both regimes. Now we need to average 
G over 77, so to define an object that may fit in the factorization formalism. 

We define (cfr eqsl3landfT0| 

f(x) = I A J dCe-^G( V ,0 = J d^e-^hGM. (32) 

I A performs the average over an 77— range of size A. If A is large enough to make the coherence region rj ~ 7 
irrelevant, the above eq 1321 is equivalent to the ordinary definition of a distribution function, with f(x) = ^(a:)! 2 - 

We assume here n = 7 /A << 1, but not by several orders. So, we expect the coherence region 7/ ~ 7 to have 
influence, although small, on the final average result. For small n it is possible to use the approximations of section 
IV. Defining 

f (x) = f(x)\ n=0 (33) 

we get (see eqfT6|) 

f(x) w f (x + in) » f ( x )+irSl = |V,(.T)| 2 +*n^^ (34) 

ax ax 

The above approximate equality derives from two calculation steps in eq J16I For the present work the interesting 
one is in the former one: 

f{x) w fo{x + in) (35) 

It requires n « 1 and x far from and 1 by an amount of size n at least. Its validity is demonstrated in [57j via 
the mapping z = log[x/{\ — x)]. However it is reasonable at intuition level: if the above conditions on n and x are 
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respected, one may guess that all the relevant things take place in a restricted area of the complex plane near the real 
axis range (n, 1 — n), and assume that the integral is dominated by singularities in this region. 

Eq l35l savs that a finite width n in S n (x — x') has the effect of introducing a finite imaginary part n in the argument 
of f (x + in). An infinitesimal imaginary part is present in / because of causality, also in no-rescattering regime: in 
the fourier transform f(x) — J d^exp{~-ix^)G(^) the lower integration limit is zero and not — oo. Since this means 
that we fourier-transform a function of the form f{x) actually means f(x + iO) (i.e. is accompanied by 

an infinitesimal damping that takes care of the potentially dangerous limit £ — > oo). A finite n makes the imaginary 
shift finite: iO — > in. It means #(£) — > 6(£)exp(— n£). 

To assume that n is small allows for the easy approximation eq !35i but the described property (the imaginary 
shift becoming finite) is likely to be more general. At intuition level, it means that a diagram where some degrees of 
freedom associated to rescattering are averaged away acquires absorption properties, and this is not surprising (see 
the discussion in the next section). 

For small n, the second step of eq|34] allows one to write the overall result as a simply readable sum, where one of 
the two terms is not affected by rescattering. As shown in |57j . the imaginary correction may lead to an observable 
asymmetry, if it affects even-odd interference terms in the partial wave decomposition of the two-point amplitude. 
This allows for the possibility to observe the effects of the discussed finite imaginary shift. 

VI. DISCUSSION. 

A. Non-hermitean hamiltonians and hidden degrees of freedom. Analogies with Feshbach's theory and the 

role of Fok-Krilov's theorem. 

An interesting point regards the analogies between the discussed argument and some known problems of nuclear 
and statistical physics. 

A nonzero T-odd observable means that an effective and not-hermitean hamiltonian is operating. "Effective" , 
because we know that strong and electromagnetic interactions do not violate T-reversal invariance. So a really 
fundamental hamiltonian cannot produce T-odd amplitudes. On the other side, if our model of the process contains 
a hamiltonian operator that does not really take into account all the relevant degrees of freedom, this hamiltonian 
becomes "effectively" non-hermitean. 

This problem has been long studied in nuclear physics (see in particular [59]). In general terms, following the 
standard Feshbach's schcmc|59], one may write a time-evolution equation for the system wavefunction tj;, and split 
the space of all the possible solutions (satisfying scattering boundary conditions) into "globally elastic" states ipE and 
"globally inelastic" tpi states. 4 

The evolution equation may be transformed into a system of two coupled equations for ipg and ipi. Using one of the 
two equations for writing ipi in terms oiipE, one may insert this relation into the other one and derive an equation for 
tpE only. The evolution of tpE is now determined by an effective hamiltonian operator that is partly non hermitean. 
The non-hermitean part is determined by "locally inelastic" processes, i.e. second order processes where an inelastic 
channel is excited as an intermediate state, to finally de-excite into the elastic channel. The Green function for the 
"reduced" problem (the problem for the globally elastic channel only) has a pole with a finite imaginary part for the 
self-energy correction: 1/(E - E' + i0) -> 1/(E - E' + iT). 

As well known in nuclear physics, these arguments give a solid formal framework for effective, non-hermitean, 
hamiltonian operators (optical potentials, etc), but are rarely useful for their direct calculation from first principles. 
However, they demonstrate that hiding degrees of freedom leads to an imaginary part in the energy-shell-projected 
amplitude. 

The factorization formalism in hard scattering imposes an a;— conservation shell for the propagating quark hole. In 
a situation where x is not conserved at all (i.e. in presence of rescattering) this creates the same situation as in the 
optical potential theory: a projection on a subspace of x— conserving states in one case, of energy-conserving states 
in the other case. 

In the hard scattering case the hamiltonian is substituted by a light-cone operator, but the general arguments could 
be repeated step by step, just after changing the name of the variables. 

The physical origin of the imaginary component of the potential, and consequently of the finite imaginary con- 
tribution to the self-energy, is interesting. It derives from the presence of several incoherent intermediate inelastic 
channels. If we consider a wavepacket with x— spread Sx, all the single plane wave components of this wavepacket are 



4 The word "global" means that we only require E(l) = E{2), not caring what happens in intermediate states. 
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elastically scattered, but coherence between them is lost, and the wavepacket shape decays in a "time" 1/Sx. This 
can be considered an application of the Fok-Krilov theorem[60], stating that a wavepacket whose frequency spectrum 
is made continuous by an interaction for t > describes an irreversible process, and has a decaying self-correlation 
function (see [4l[ for extensive discussion about this). 



B. Partonization/Decorrelation process within known models 



In the example of section V we have considered an ideal limit for the chaotic hadron-to-parton evolution, that could 
be named "ideal partonization process" , where the probabilistic features of the initial quark state with respect to x 
are perfectly conserved, although phase features are lost. 

In terms of known quark-hadron diagrams, an example of a quark state in a hadronic phase is the hadron-quark- 
spectator vertex employed in the spectator model by [61| . Such a splitting diagram describes a quantum fluctuation 
of a hadron into a quark+spectator pair. This fluctuation is perfectly reversible to the original hadron state. 

The presence of a soft cutoff A ~ Rhadron i n the vertex form factor guarantees that no momentum component k > 
A crosses the vertex, so that correlation is present up to a distance Sz ~ R. 

When the splitting takes place under the action of initial state interaction with another hadron, the process is 
incoherent and presents a certain degree of irreversibility. Then we may speak of "partonic phase" . 

A typical example of such a process is a diagram for the calculation of the Sivers function, where a gluon exc hang e 
is added to the previous splitting process (it may be found in several of the above quoted models, starting from [2l|). 
The added gluon loop involves momenta up to the hard scale, and this destroys soft states. 

At the lowest order, an example of transition between the two situations is a cut diagram where the previous two 
examples interfere. This interference graph is the basic ingredient of several calculations of the Sivers function, and 
literally reproduces the light-cone evolution from one of the two phases to the other one (taking the causality factor 
£?(£) into account). 



VII. CONCLUSIONS 



In this work an example has been shown for the correlator of a hadron state, that presents time-evolution properties. 
It initially consists of a pure state, and finally of a completely incoherent probabilistic mixture. The a;— distribution 
is the same in both cases. This evolution is caused by rescattering, and is visible as far as the additional degrees 
of freedom associated with rescattering are explicitly visible. After these degrees of freedom have been integrated, 
the remaining correlator shows modified analytical properties, consisting in a finite imaginary shift of the main cut 
singularity of the process. 
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